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Motivated by the experiments on DNA under torsion, we consider the 
problem of pulling a polymer out of a potential well by a force applied to one 
of its ends. If the force is less than a critical value, then the process is activated 
and has an activation energy proportinal to the length of the chain. Above 
this critical value, the process is barrierless and will occur spontaneously. We 
use the Rouse model for the description of the dynamics of the peeling out and 
study the average behaviour of the chain, by replacing the random noise by 
its mean. The resultant mean-field equation is a nonlinear diffusion equation 
and hence rather difficult to analyze. We use physical arguments to convert 
this in to a moving boundary value problem, which can then be solved exactly. 
The result is that the time tpo required to pull out a polymer of segments 
scales like A^. For models other than the Rouse, we argue that tpo A^^*^ 

I. INTRODUCTION 

Considerable attention has been paid to the replication of the DNA molecule |1| . The first 
step in the replication is the unzipping of the two strands. Usually, this is caused by enzymes 
and mechanical force at the molecular level is involved in the action of these enzymes. 
Recently, micromanipulation techniques have been developed to study single molecules of 
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DNA under stress 0. Most such single molecule experiments have been on the stretching 
properties of DNA, but very recently, the response of the DNA to external torques has also 
been studied 0,^. Mechanical unzipping of the two strands of the DNA, in the absence 
of enzymes has been carried out p. There have been interesting investigations in to the 
theory of stretching the DNA |^ and its denaturation under a torque (see Fig. |I]) In 
an interesting recent e-print, Bhattacharjee has suggested a minimal model to study the 
unzipping. He treats the DNA as consisting of two flexible interacting elastic strings which 
are bound together by an attractive interaction. The two are tied together at one end and at 
the other end, there are forces acting on the two strands, trying to separate them. Winding 
is ignored in his analysis. He assumes that in relative coordinates, the Hamiltonian for the 
unzipping problem is equivalent to that of for a single chain, subject to an attractive short 
range potential well with a force acting at one of its ends. It is shown that the polymer 
can be pulled out only if the force exceeds a critical value. The model is of interest in other 
contexts too. For example, one can imagine constructing micro-fabricated channels having 
different depths in different regions to trap the polymer by entropic barriers (see the recent 
experiments by Han et. al. [Q) and then pulling it out from the traps by the application of 
a force at one end. Another possible application is the desorption of a long chain molecule 
adsorbed on a surface by pulling at one end. Our aim in this paper is to investigate the 
dynamics of the escape of a chain from a potential well as a result of applying a force at 
one of its ends. We consider only the one dimensional problem and our approach is much 
less sophisticated than that of 0, as we do not take in to account the excluded volume 
interaction. We bring out the physical reason for the existence of the critical force. For 
describing the dynamics, we use the simple Rouse model. We find that tpo, the time required 
to pull out a chain molecule of N segments scales as A^^. 

In the section [I^ we give an outline of the statistical mechanics of a polymer, trapped 
by an attractive potential, with a force acting at one end. Section |IT1| outlines the Rouse 
model that is applicable for a polymer subject to a force, and analyze why the force has to 



exceed a critical force for pulling it out. In section we argue that the problem can be 



approximated by a moving boundary value problem and then solve the problem. Finally, 
section |V] summarizes our conclusions. 



II. STATISTICAL MECHANICS OF A TRAPPED POLYMER WITH A FORCE 

AT ONE OF ITS ENDS 

In this section we consider the statistical mechanics of a polymer in one dimension, 
trapped in an attractive potential V{R), that is located near R = (see |^ for the nature of 
the potential) and subject to a force F at one of its ends. We take the number of segments 
in the polymer to be N. The partition function for a free polymer, that is not subject to 
any force, but having one end at i? = may be calculated as Z = / dRG{R, 0, A^) where 
G{R, 0, N) is the propagator for the unforced molecule, which obeys the differential equation 
( - f & + k^ViR))GiR, Ro, N) = SiN)SiR - Ro) (see |g). lis the Kuhn length for 
the chain and R and Ro denote the positions of the two ends of the chain. 

In the presence of a force F acting at the end at R, the partition function becomes 

RF 

dRG{R,0,N)e''BT (1) 
Using the ground state dominance approximation for the propagator, valid for long polymers 



10| , we can write the above as: 

poo rf-EqN 

dRro{R)MQ)e^^ (2) 



where ipo{R) is the lowest eigen-function of the operator H = — f-^2 + having the 



eigenvalue EQ/{kBT). Now, for large values of R, the eigen-function ipo{R) ~ exp{—^p^^R) 



~ exp(— y p^^'y -R), where in the last step we have neglected the "zero point energy" contri- 
bution to the lowest eigenvalue and approximated it by — Vq, where Vq is the depth of the 



potential well. On using this in the equation (0), we realize that for F > Fc = y^^^r-^, the 
contribution from large values of R make the integral diverge, indicating that for F > Fc, the 
polymer gets pulled out of the hole. For F < F^, the partition function is finite, indicating 
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that the molecule remains trapped. In this analysis, we have neglected the zero-point energy, 
which in the polymer problem, is equivalent to neglecting the effect of thermal fluctuations. 
Thermal fluctuations make Fc slightly lower 



III. THE ROUSE MODEL AND THE FREE ENERGY PROFILE FOR THE 

PULLED POLYMER 

The process that we study is shown in the Fig. ^. For its description,we use the continuum 
limit of the Rouse model, discussed in detail by Doi and Edwards fl^, which has recently 



been used to study the related problem of escape of a polymer over a barrier |TTI]. The 



chain is approximated as a string, with segments (beads) labelled by their position n along 
the chain, n is taken to be a continuous variable, having values ranging from — to 
(this is convenient for the present problem). The position of the n*^ segment in space (one 
dimension) is denoted by i?(n, t), where t is time. In the Rouse model, this position undergoes 
overdamped Brownian motion and its time development is described by the equation 

C— ^ = - ^'(^K t)) + /(n, t)- (3) 

In the above, C is a friction coefficient for the n^'^ segment. The term an"'*'' comes from 
the fact that stretching the chain can lower its entropy and hence increase its free energy. 
Consequently, the parameter m = SksT/l'^ (see Doi and Edwards |TD[, Eq. (4.5). They use 



the symbol k for the quantity that we call m) . V{R) is the free energy of a segment of chain, 
located at the position R and V'{R) = dV{R)/dR. f{n,t) are random forces acting on the 
^th segment and have the correlation function (/(n, t)/(ni, ti)) = 2(kBT6{n — ni)6(t — ti) 
(see Eq. (4.12)). The deterministic part of the Eq. (|^), which will play a key role in 



our analysis, is obtained by replacing the random noise term in (^), by its mean. As we 
are applying a force at one end of the polymer (at n = 0) and as the other end is free, the 
boundary conditions to be satisfied are 
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'dR(n,t) 



[ dn 

The above equations may equivalently be written as 



= 0. (5) 

n=-N 



where we have defined the free energy functional 

E[R] = dn I ^ + V{R) - FR{n)5{n + e) \ (7) 



e (— ^ 0), is a small positive number. The process of pulling out may be made physically 
clearer by analyzing this free energy functional. For this, we consider the free energy hyper- 
surface that results from the Eqs. (^) (j^) and (^. We find the extrema on this surface by 
putting j^E[R\ = 0, which leads to 

m—r = V'iR) for n< (8) 

This has to be solved, subject to the conditions of the Eqs. (Q) and (^. Imagining and 
speaking of n as time for the rest of this section, we realize that this is just Newton's 
equation for a particle of mass m, moving in the upside down potential —V{R). In the Fig. 
^, we show this potential. The boundary condition of Eq. (H) means that the particle has to 
start at the time n = —N, with zero velocity and end at the time n = 0, with the velocity 
F/m. This means that at the end of its path, its kinetic energy is At n = —N, the 
particle can at the most start at the top of the potential, where its total energy is Vq. If it 
obeyed the equation (|]), then at n = 0, its kinetic energy has to be less than this amount. 
Thus we get the condition that Vo > ^ (or equivalently, Fc > F) for the free energy surface 
to have a local minimum). If this condition is satisfied, then for any finite value of A^, we 
can find two solutions to the equation (H) satisfying the boundary conditions (this is true for 
potentials of the shape shown in the figure). These two solutions are shown in the Fig. ^ In 
the first solution (shown as line AB in the figure), the particle starts near the maximum of 
—V{R), at the time n = —N and takes a very long time to move away from the maximum. 



and eventually arrives at the point B at the final time n = with exactly the right kinetic 
energy of As this corresponds to a configuration for the polymer where almost all its 
(assumed to be large) segments are near the minimum of the potential well V{R), we 
can estimate the free energy for this configuration to be ~ — A^Vq. The second solution, also 
shown in the Fig. ^, CD, corresponds to the particle starting from rest from an appropriate 
value of R (point C in the Fig. ^ at the time n = —N and arriving at n = with the 
kinetic energy The particle spends only relatively small amount of time in the region 
where the potential energy —V{R) is non-zero. One easily estimates the free energy of the 
polymer to be ~ Thus, for a polymer of A^ segments, there is a barrier of height 

A^(Vo — ^) = N(F^ — F'^)/2m, to be overcome for the pulling the polymer out. As F ^ 
this barrier is reduced and if F > Fc then the free energy surface has no local minimum - 
its minimum is at infinity - this means that the pulling out is barrierless and would happen 
spontaneously. A pictorial representation of this scenario is given in the Fig. ^. 

IV. THE DYNAMICS 

We now consider the case with F > F^ so that the pulling out of the chain has no 
activation barrier. We are interested in only the average behavior of the segments, in the 
process of pulling out. This may be obtained by analyzing the deterministic equation 

^dR(n,t) d'^R(n,t) 

where we have replaced the random term f{n,t) by its average value, viz., 0. We now wish 
to analyze this equation subject to the two boundary conditions of the Eqs. (|D and (H). 
This equation, is just a nonlinear diffusion equation. 

The equation (|^) is very difficult to solve for any realistic V{R). However, we can use 
the following physical picture. We take V^(-R) to be a short ranged well having a depth of 
Vq, and to have a rather sharp boundary like the one in the figure ^ The change over from 
—Vq to the fiat region where V{R) = occurs over a region of width w, which is assumed to 
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be small. This means that V'{R{n,t)) resembles a Dirac delta function, which implies that 
^^q"'^^ has to change very rapidly, across this region, and would resemble a step function. If 
we consider the process of escape from the well, then at any time, the units of the polymer 
with n < riiit) would have R{n, t) = 0, and the segments from riiit) to rioit) would be spread 
across the region of width w. We shall refer to no{t) as the peeling point. On multiplying 
the equation (||) by ^^q^'^^ and integrating from to no{t), we get 



m fdR{n,t)Y rMt) dR{n,t) dR{n,t) 



/ nM) •^"^ 



dn I J riiit) dt dn 



dn (10) 



In arriving at the Eq. (|ToD we have used the fact that [ ^q^'^^ j = and V {R{ni{t))) 



— Vq. If the width w is small, then we are justified in neglecting the right hand side of the 
Eq. ([lOD and then we find that ( ^^g"'*'* ) = ^! which on taking the square root may be 
written as 

( dR{n,t) \ ^ ^ ^^^^ 
where Fc = ^/2mVo• (If we picture the process as in Fig. ^, then we need not worry about 



the case where Fc = — y^mVo). Now, all that we have to do is to solve the diffusion equation 
(y), outside the well region, subject to the above boundary condition. Outside the well, the 
equation (|^) becomes the simpler diffusion equation 

, dR{n,t) d^R{n,t) 

which is a linear equation. As the range of the potential well is small in comparison with 
the total length of the polymer, we can take 

R{n,t) = for n < no{t) (13) 

Thus we now have a moving boundary value problem: we have to solve the Eq. (|1^) for 
R{n,t), and determine Uoit), such that the conditions (|), ( [Tl| ) and ( [T3| ) are satisfied (see 
Fig. ^ for a pictorial representation of the process). In thus simplifying the equation, our 
assumption of neglecting the right hand side of the Eq. (|lOl) is crucial. Physically, what 



is being done is that near the point of breaking away from the potential, the segments are 
assumed to be at equilibrium - that means the time scale for peeling the polymer away from 
the well is slow in comparison with the time scale for establishment of local equilibrium near 
the peeling point. 

Usually, moving boundary value problems are difficult to solve. However, in this case, it 
is possible to find similarity solutions to these equations. For this, we put R{n,t) = \/tp{z), 



where z — n/{a^/t), with a — sj^- 



(Pp{z) dp{z) 



dz'' 



+ Z' 



dz 



-p{z)^0 



(14) 



Note that z e (—00, 0). The boundary conditions become 



' dp{z) 
dz 



z=0 



aF 
m 



(15) 



' dp{z) 
dz 



aF, 



m 



(16) 



and 



p{zo) = 



(17) 



where Zo is related to no{t) by 



oit) = ZoaVt 



These equations can be solved to find p{z) and Zo- The result is: 



a 



Fc -F)\\l + ^ {{Fc - F) Erf (^) + F Erf (^) 



mErf(^) 

where Zo{< 0) is the solution of the equation 

- Erf(^)e . .0 = 

Note that Zq is a function of ^'^ p'^"-' only. Using all the above results, we get: 



(18) 



(19) 



^.Kl^^n.^. il^ , &-fi^ + „ >^Erf(^„ ./X) I (20) 
F V 2t V2mt FErf(^)\^V7r V 2mt ^^2 \ 2mt' ' ^ ' 

To demonstrate the nature of this solution, we note that the right hand side is a function of 
n \J 2^ alone and make a plot of the left hand side against this variable, in the Fig. |^. The 
important conclusion that we can draw from the above is that the point of detachment of 
the polymer from the well is given by the equation no{t) = zq^J^^, where zq is the solution 
of the equation (|l^) and is t independent. As the left hand side of the equation ([19| ) is 



positive, one has a solution for zq only if ^^^^"^ > and again this just means the existence 
of the critical force (a similar conclusion is arrived at by Bhattacharjee 0). His answer is 
slightly different, as his analysis takes the noise term too in to account). The pulling out of 
the polymer is complete when no{t) = —N. Hence the time required to pull out the chain is 
given by 

tr. 



N ^ 2 



Zoa 



For forces slightly exceeding the critical force, one can solve the Eq. ([T9|) and get Zq = 

1 /2 

— ■ Hence we conclude: tpo ~ A^^. Further, for small (F — Fc), 

(21) 



V. CONCLUSIONS 

We now summarize our results and point out the drawbacks of our analysis. Using the 
Rouse model, we find that for a polymer trapped in a potential well, subject to a force at one 
of its ends, there exists a critical force F^ and below this critical force, the pulling out is an 
activated process, with an activation energy proportional to the number of segments in the 
chain. So the pulling out would not happen. For forces higher than F^, there is no barrier 
and the process would occur spontaneously. The dynamics of the pulling out is described 
by a non-linear Rouse equation, which is rather difficult to solve. We find an approximate 
solution by (a) confining ourselves to the average behavior of the positions of the segments 
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and (b) introducing a local equilibrium assumption, which enables us to convert the non- 
linear equation to a moving boundary value problem. This moving boundary value problem 
has a similarity solution and this leads to the conclusion that the time required to pull the 
chain out scales like iV^. 

Obviously, there is scope for improvement of our analysis. One would like to perform 
a more detailed analysis, including the effects of noise. Also, the effects of hydrodynamic 
interactions and the excluded volume interactions have to be included. Inclusion of these 
may be done by modifying the equation || to 

^aR^ = m / dnMn - ^^)^^^^ " V'i^i^^ ^)) + ^) 
(see the equation (4.50) of Ref. fl^). In general h{n) ~ n"" with f < 1 (different values can 
be used for u to account for Zimm dynamics or excluded volume interactions - see equations 
(4.49) and (4.69) of Ref. [0). With the boundary conditions same as in our equations (H), 
(0) and (|13D, then one finds that tpo ~ iV^+". 

Our results above, however, are only the first few steps in the analysis of this very 
interesting problem. The most important next step, obviously, is to go beyond our mean 
field description by including the fluctuations. It seems likely that fluctuations can modify 
the long time behavior qualitatively. 
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FIG. 1. The forced unzipping of DNA. 
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FIG. 3. The upside down potential —V{R). The hnes AB and CD represent the solutions of 
the Eq. ^ corresponding to the initial and the transition states respectively. The corresponding 
configurations of the polymer are represented by the coiled curves above them. TS stand for 
transition state. The kinetic energy at the n = end is ^ and this must be less than the height 
of the potential Vq, if the extrema are to exist. 
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FIG. 4. The free energy hypersurface. For F < Fc it has both a local maximum and a minimum, 

n(f'^-f^) 

while for F > Fc no local extrema exist. Eact is the activation energy and is equal to ^ ^ — 
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FIG. 5. The moving boundary value problem. The force F is applied at n = 0. n varies from 
—N to 0. At a given time t, the portion with > n > no{t) has been pulled out from the well. 
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FIG. 6. Plots of J^^^ against 
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